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We investigated the diversity and community composition of denitrifying bacteria in surface water from the Yellow River 
estuary. Our results indicated that the diversity of the denitrifying community in freshwater based on the nirK gene was higher 
than that in seawater. Furthermore, phylogenetic analysis suggested that the bacteria community could be distributed into eight 
clusters (Clusters I to VIII). Redundancy analysis (RDA) revealed that community compositions were related to multiple 
environment factors, such as salinity and nitrate concentration. The results of the present study have provided a novel insight 
into the denitrifying community in water columns in estuaries. 
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MICROBES AND 
ENVIRONMENTS 



Nutrient removal has played an important role in pre- 
venting the eutrophication of receiving waters (22). 
Denitrification, as an effective way to remove nitrogen, is an 
alternative anaerobic respiration process that removes nitro- 
gen via the microbial stepwise reduction of N0 3 " to gaseous 
products (NO, N 2 0, N 2 ) (5). Nitrite reductase is the rate- 
limiting enzyme in denitrification, and catalyzes the step from 
nitrite reduction to nitric oxide. Thus, it has commonly been 
used as a molecular marker of denitrifying bacteria. Two 
main classes of functionally-equivalent nitrite reductase have 
been identified in denitrifying bacteria: a copper-containing 
NirK enzyme and cytochrome cd l NirS nitrite reductase (2, 
16, 25). The nirK gene, which encodes nitrite reductase, has 
extensively been used in recent years to clarify the composi- 
tion of the denitrifier community and diversity in the water 
columns of freshwater and seawater (10, 14, 17, 21). 

An estuary is a complex ecosystem that receives extensive 
river discharges of various terrestrial and anthropogenic 
materials, such as nutrients and pollutions (9). Terrestrial 
inputs and physiochemical parameters have been shown to 
have a marked impact on the diversity and community 
composition of denitrifying bacteria in the water column of 
an estuary. The Yellow River estuary is located in the eastern 
coastal area of China, and high concentrations of nitrogen 
are a feature in this water column. The importation of 
nitrogen has recently been shown to be increasing in the 
Yellow River estuary (13). However, the diversity and 
community composition of denitrifying bacteria in surface 
water in the Yellow River estuary remain unknown. In this 
study, we investigated the diversity and distribution of 
denitrifying bacteria in the water column in the Yellow River 
estuary and revealed relationships between the denitrifying 
community and environmental parameters. 

Four samples were collected from the Yellow River 
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estuary on October 21, 2010 (Fig. SI). Samples from sites A 
and B were taken from freshwater sources, whereas samples 
from sites C and D were from seawater sources. Water sam- 
ples at each site were taken in triplicate at a depth of 0.5 m by 
a water sampler (Wildlife Supply Company, USA). A 1 L 
water sample from each site was filtered through 0.22 |um 
millipore filters. These filters were stored at -80°C until DNA 
was extracted. The physicochemical variables at each sam- 
pling site were measured three times and the average values 
are shown in Table 1 . Significant differences were observed 
in the salinity in all sampling sites, ranging from 2.3 to 26.7 g 
L" 1 . The concentrations of total nitrogen, dissolved oxygen, 
and nitrate were markedly higher in freshwater than in 
seawater samples. In contrast, the concentrations of chemical 
oxygen demand were lower in freshwater than in seawater 
samples. 

The genomic DNA of each sample was extracted using 
the E.N.Z.A.™ Water DNA Kit (Omega, USA) according to 
the manufacturer's instructions. Fragments of the nirK and 
nirS genes were amplified using the primer pairs FlaCu- 
R3Cu for nirK (6) and nirSlF-nirS6R for nirS (3). The PCR 
amplification conditions are shown in Table SI. No nirS PCR 
products were obtained in any of the four samples by repeated 
PCR. The purified PCR products were ligated into the 
pMD18-T simple vector (TaKaRa, Japan), and then trans- 
formed into Escherichia coli DH5a competent cells to con- 
struct the gene libraries. Approximately 80 colonies were 
selected from the clone library of each sample. The clones in 
each library were screened by restriction fragment length 
polymorphism (RFLP). The PCR products (8 |uL) were added 
to 20 uL reactions and incubated at 37°C for 2 h, containing 
1 U each of the enzymes Haelll and Msp\ (TaKaRa, Japan) 
and 2 |uL buffer (24). The representative clones were then 
selected for sequencing in Majorbio Biomedical Technology 
Co. Ltd. (Shanghai, China). 

Amino acid sequences sharing 95% similarity were clus- 
tered into a single operational taxonomic unit (OTU 0 95 ) by 
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Table 1. Physicochemical properties of water column samples collected from the Yellow River estuary 



Sample 


Tem a (°C) 


pH 


SalMgL" 1 ) 


TN C (mg L 1 ) 


N0 3 "-N d (mg L 1 ) 


TP e (mg L 1 ) 


DO f (mg L 1 ) 


COD§ ( mg L i) 


A 


13.9 


8.38 


2.3 


11.0* 


2.36* 


0.072 


9.5 


20.2 


B 


13.9 


8.45 


3.0 


5.7 


3.35* 


0.115 


9.5 


18.2 


C 


14.1 


8.14 


24.7* 


3.1 


0.11 


0.079 


7.6 


63.0* 


D 


14.1 


8.10 


26.7* 


4.5 


0.04 


0.079 


7.5 


53.6* 



a, temperature; b, salinity; c, total nitrogen; d, nitrate; e, total phosphorous; f, dissolved oxygen; g, chemical oxygen demand; *, significant differ- 
ences (P<0.05) according to Duncan's multiple range test using SAS version 9.1 software. 



Table 2. Community diversity and predicted richness of NirK sequences from each sampling site 



Sampling site 


No. of clones 


No. ofOTUs a 


C(%) b 


^ACE C 


c d 
^Chaol 


D e 


H r 




A 


63 


24 


77.8 


59.0 


55.1 


0.117 


2.56 


0.80 


B 


60 


23 


78.3 


46.3 


49.0 


0.067 


2.79 


0.89 


C 


50 


14 


88.0 


20.1 


19.8 


0.137 


2.19 


0.83 


D 


56 


14 


92.7 


19.0 


17.4 


0.116 


2.27 


0.86 



a, Number of OTUs defined by the furthest neighbor algorithm in MOTHUR at 95% similarity; b, Coverage=l-(/7/A0 x 100, where n t refers to the 
number of clones appearing only once in each library and TV refers to the total number of clones in the library; c, Richness-based coverage estimator; 
d, the estimated richness; e, Simpson's index; f, Shannon- Wiener index; g, Evenness index. 



MOTHUR software (18). Phylogenetic trees were con- 
structed by the MEGA 5.1 program (19) using the neigh- 
bor-joining method and maximum composite likelihood 
model. The diversity indexes were calculated by Biodap 
software (20). Redundancy analysis (RDA) was performed in 
CANOCO 4.5 for Windows (1). 

In total, 229 clones of nirK genes were analyzed, with 
63, 60, 50 and 56 clones being obtained from sites A, B, C, 
and D, respectively. The numbers of OTUs in each library 
were 24, 23, 14, and 14, respectively, and the coverage of 
each library varied from 77.8% to 92.7% (Table 2). The 
library of site A had the highest richness based on ACE and 
Chaol, while the Shannon- Weiner index and Simpson's 
index indicated that the diversity of site B was higher than 
that of other samples. We also found that the rarefaction 
curves for seawater samples were markedly flatter than those 
for freshwater samples (Fig. S2). These results revealed that 
diversity and richness were markedly higher in freshwater 
than in seawater samples. 

The NirK compositions and relative ratios differed among 
different sampling sites (Fig. 1). OTU30 and OTU36 were 
the dominant OTUs in sites A and B. Nevertheless, most 
dominant clones in sites C and D were significantly different 
from those in sites A and B. The most dominant OTU was 
OTU16 in site C, followed by OTU17, whereas OTU13 and 
OTU6 were the dominant OTUs in site D. Detailed data are 
shown in Table S2. No common OTU was shared in the four 
sampling sites (Fig. 1 and Fig. S3). Community composition 
analysis and Venn diagrams revealed that OTU5, OTU30, 
and OTU36 appeared in sites A, B, and C, but were absent 
in site D. Furthermore, OTU27 existed in all sampling 
sites, except for site C. Furthermore, twelve OTUs were 
shared in more than one site, and 47 OTUs were detected in 
only one site (Table S2 and Table S3), which indicated the 
representive OTUs of each site. 

A neighbor-joining (NJ) phylogenetic tree based on amino 
acid sequences was generated from the 59 NirK OTUs in all 
sampling sites (Fig. 2). The NirK sequences from these 
samples were grouped into eight Clusters (I to VIII). Cluster I 
contained sequences from four sampling sites, the reference se- 
quences of which came from various environments including 
lake water, sewage, and sediment (8, 11, 12). Thus, these 
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Fig. 1. Relative abundance of OTUs in sampling sites A, B, C, and D. 

results indicated that Cluster I may be a ubiquitous denitrify- 
ing group. Clusters IV and VI also contained sequences from 
all sampling sites (mainly from seawater samples). These 
sequences were similar to those in two lakes (lakes PluBsee 
and Schohsee), the Baltic Sea (10), and the San Francisco 
Bay estuary (12). Cluster III contained sequences from all 
sampling sites except site D, which was closely related to 
clones previously described from activated sludge (7, 15). 
The main sequences of Clusters II and V were from fresh- 
water sources (sites A and B), and were closely related to 
those from Lake Kinneret (8), activated sludge (7, 15), and 
the San Francisco Bay estuary (12). The NirK amino acid 
sequences deduced from nirK gene sequences in the Yellow 
River estuary were closely related to sequences from differ- 
ent habitats, which indicated that these clones were widely 
distributed denitrifying bacteria and may be able to flexibly 
adapt to various environments. 

RDA analysis was employed to determine the influence of 
environmental factors on the mriC-harboring denitrifier com- 
munity (Fig. 3). The first and second dimensions explained 
82.9% and 13.4% of the total variance, respectively. RDA 
analysis revealed that the community compositions of denitri- 
fying bacteria were related to multiple environment factors, 
such as salinity and nitrate concentration. The results 



nirK-Harboring Bacteria in the Yellow River Estuary 
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OTU37 (1 A clone, 2 B clones) 
Lake Kinneret water (ABW75 148) 
OTU36 (16 A clones, 9 B clones, 2 C clones) 
OTU38 (1 B clone) 
r OTU39 (1 A clone) 
L Norway sewaee (ABZ86195) 
OTU50 (1 B clone) 
Agricultural plots soil (ABI 19341) 
OTU29(l A clone) 

■Bradyrhizobium sp.STM 3809 (ZP 09427849) 

- OTU28 (1 B clone) 

OTU27 (4 A clones, 5 B clones,4 D clones) 
San Francisco Bay estuary sediment (ADM94059) 
Sweden activated sludge (ABA54053) 
OTU34 (4 B clones) 
OTU33 (I A clone,4B clones) 

■ OTU35 (1 A clone) 

^Tl— Two lakes and Baltic Sea water (ABR291 10) 
OTU31 (1 A clone) 
OTU32 (1 A clone) 
•Paracoccus mR-26901 (CAJ76774) 
OTU21 (2 A clones,3 B clones) 
OTU30 (14 A clones, 10 B clones, 1 C clone) 
.„ Activated sludge (BAE43773) 
Sweden activated sludge (ABA54018) 

- OTU57 (1 B clone) 
OTU56(l A clone) 
OTU58 (1 A clone) 
-OTU52(l A clone) 

Citromicrobium bathyomarinumJL354 (ZP 06863130) 

- OTU51 (2 B clones, 3 C clones) 
OTU55 (1 A clone) 

ioo, OTU53 (2 A clones) 

Activated sludge (BAD36893) 
■Agrobacterium tumefaciens (ZP 12909478) 
OTU54 (1 A clone) 

■ Activated sludge (BAE43799) 
^.Pseudomonas mendocina NK-01 (YP 004377929) 
" ' OTU40 (1A clone,8 D clones) 

Pseudomonas mendocina (AAL82514) 
OTU41 (1 A clone) 

- OTU42 (2 A clones) 
OTU43 (1 B clone) 

- Sweden arable soil (ABC02745) 
OTU48 (5 C clones) 

- San Francisco Bay estuary sediment (ADM941 15) 
OTU49 (4 C clones) 

- OTU47 (5 D clones) 
Silicibacter sp. TrichCH4B (ZP_05742967) 

OTU46(l B clone) 
OTU45 (1 C clone) 

■ OTU44 (1 A clone) 

Two lakes and Baltic Sea water (ABD36616) 
OTU25 (1 B clone) 

■ Lake Kinneret water (CAJ44502) 
OTU23 (3 B clones) 

OTU24 (1 B clone) 

San Francisco Bay estuary sediment (ADM94042) 
OTU26 (2 B clones) 
OTU20 (4 B clones) 
-OTU22(l B clone) 

- Mexico soil (AER93073) 
■ OTUl 9 (4 A clones, 1 B clone) 
1 San Francisco Bay estuary sediment (ADM94044) 

OTU16 (13 C clones) 
OTU17 (12 C clones) 

Two lakes and Baltic Sea water (ABD36591) 
OTUl 8 Q C clones) 

- OTU15 (2 C clones) 

- Arabian Sea water (ABW35183) 
OTU7 (1 B clone) 
OTU6 (10 D clones) 

San Francisco Bay estuary sediment (ADM93845) 
Two lakes and Baltic Sea water (ABD36652) 
OTU5 (3 A clones, 1 B clone, 1 C clone) 
• OTU9 (1 D clone) 
OTU8 (1 D clone) 

San Francisco Bay estuary sediment (ADM93918) 

- OTU14 (1 C clone, 1 D clone) 
Two lakes and Baltic Sea water (ABD3661 1) 

- OTU13 (1 C clone, 13 D clones) 
San Francisco Bay estuary sediment (ADM94192) 
OTU12(l D clone) 

■ OTU10 (2 D clones) 
OTUl 1 (1 D clone) 
Jiulong River sediment (ADR31859) 

-OTU4(l A clone) 

- San Francisco Bay estuary sediment (ADM93852) 
Shewanella loihicaYMA (YP 00 1095460) 





Cluster I 



Cluster II 



Cluster III 



Cluster IV 



Cluster V 



Cluster VI 



OTUl (2 D clones ) 
- OTU2 (3 D clones) 
OTU3 (4 D clones) 

r OTU59 (1 C clone) 

loo 1 - San Francisco Bay estuary sediment (ADM93913) 
Moraxella catarrhalisRHA (YP_003626300) 



Cluster VII 



Cluster VIII 



0.1 



Fig. 2. Phylogenetic analysis of denitrifying bacteria based 
on NirK sequences obtained in the Yellow River estuary. A 5% 
cut-off in the amino acid sequence was used to define OTUs by 
MOTHUR. The neighbor-joining method was used and boot- 
strap analysis was performed with 1,000 replications. Bootstrap 
values above 50 were indicated at branch points. The aniA gene 
from Moraxella catarrhalis (YP 003626300) was used as an 
outgroup. 



obtained also indicated that Clusters I, II, and V positively 
correlated with total nitrogen, nitrate, dissolved oxygen, total 
phosphorus, and pH, and negatively correlated with tem- 
perature, salinity, and chemical oxygen demand. The oppo- 
site results were obtained for Clusters IV and VI (Fig. 3 and 
Table S4). The relationship among sampling sites, commu- 
nity composition, and environmental factors indicated that 
samples from sites A and B were more similar than those 
from sites C and D. Dang et ah (4) reported a relationship 
between denitrifying bacteria and the surrounding environ- 
ment in the Jiaozhou Bay, and the different environmental 
adaptation strategies of various denitrifying bacteria was 
subsequently proposed. We demonstrated in the present study 
that salinity significantly influenced the diversity and com- 
munity composition of denitrifying bacteria. On one hand, 
we found that the diversity of the denitrifying community 
inversely correlated with salinity (Table 1 and Table 2). 
Similar results have also been reported in previous studies on 
different habitats (17, 23). On the other hand, the community 



compositions of denitrifiers were distinct due to their salini- 
ties (Fig. 2 and Fig. 3). This result may have been caused by 
the selection effect of salinity. In a previous study, denitrify- 
ing communities in the freshwater of Lake Kinneret were 
shown to differ from those of marine habitats, suggesting the 
differentiation of marine and freshwater denitrifying bacteria 
(8). These results revealed that environmental factors may 
play critical roles in controlling the diversity and spatial 
distribution of the denitrifying community. 

In summary, this study showed that denitrifying bacteria 
existed in surface water and were very diverse in the Yellow 
River estuary. The community compositions and relative 
ratios of bacterial communities markedly changed with the 
different sampling sites. To the best of our knowledge, this is 
the first study to examine the denitrifying bacteria community 
in the Yellow River estuary, and the results obtained have 
provided a novel insight into the denitrifying community in 
the estuary, especially in surface water. 

The representative sequences of nirK fragments reported in 
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Fig. 3. The RDA ordination plot for the relationship between the 
distribution of Clusters and environmental factors in the Yellow River 
estuary. Correlations between clusters or environmental factors and 
RDA axes are represented by the length and angle of arrows. A, B, C, 
and D represent the four sampling sites. 

this study have been deposited in GenBank under the acces- 
sion numbers KF 143 898 to KF 144045. 

Acknowledgements 

This study was financially supported by grants from the National 
Natural Science Foundation of China (No. 41306150), Promotive 
Research Fund for Excellent Young and Middle-Aged Scientists of 
Shandong Province, China (No. BS2012HZ01 1), A Project of 
Shandong Province Higher Educational Science and Technology 
Program, China (No. J10LC09), and Open Funding Project of Key 
Laboratory of Marine Biogenetic Resources, SOA, China (No. 
HY201205). 

References 

1. Braak, C.J.F. ter., and P. Smilauer. 2002. CANOCO Reference 
Manual and CanoDraw for Windows User's Guide: Software for 
Canonical Community Ordination (version 4.5). Microcomputer 
Power, Ithaca. 

2. Braker, G., J.Z., Zhou, L.Y. Wu, A.H. Devol, and J.M. Tiedje. 2000. 
Nitrite reductase genes (nirK and nirS) as functional markers to 
investigate diversity of denitrifying bacteria in Pacific northwest 
marine sediment communities. Appl. Environ. Microbiol. 66:2096- 
2104. 

3. Brettar, I., E.R. Moore, and M.G. Hofle. 2001. Phylogeny and 
abundance of novel denitrifying bacteria isolated from the water 
column of the central Baltic Sea. Microb. Ecol. 42:295-305. 

4. Dang, H., C. Wang, J. Li, T. Li, F. Tian, W. Jin, Y. Ding, and Z. 
Zhang. 2009. Diversity and distribution of sediment w'nS-encoding 
bacterial assemblages in response to environmental gradients in the 
eutrophied Jiaozhou Bay, China. Microb. Ecol. 58:161-169. 

5. Dell, E.A., D. Bowman, T. Rufty, and W. Shi. 2010. The community 
composition of soil-denitrifying bacteria from a turfgrass environ- 
ment. Res. Microbiol. 161:315-325. 

6. Hallin, S., and P.E. Lindgren. 1999. PCR detection of genes encoding 
nitrite reductase in denitrifying bacteria. Appl. Environ. Microbiol. 
65:1652-1657. 



Li et al. 



7. Hallin, S., I.N. Throback, J. Dicksved, and M. Pell. 2006. Metabolic 
profiles and genetic diversity of denitrifying communities in activated 
sludge after addition of methanol or ethanol. Appl. Environ. 
Microbiol. 72:5445-5452. 

8. Junier, P., O.S. Kim, K.P. Witzel, J.F. Imhoff, and O. Hadas. 2008. 
Habitat partitioning of denitrifying bacterial communities carrying 
nirS or nirK genes in the stratified water column of Lake Kinneret, 
Israel. Aquat. Microb. Ecol. 51:129-140. 

9. Kennish, M.J. 2002. Environmental threats and environmental future 
of estuaries. Environ. Conserv. 29:78-107. 

10. Kim, O.S., J.F. Imhoff, K.P. Witzel, and P. Junier. 2011. Distribution 
of denitrifying bacterial communities in the stratified water column 
and sediment-water interface in two freshwater lakes and the Baltic 
Sea. Aquat. Ecol. 45:99-112. 

11. Mao, Y., L.R. Bakken, L. Zhao, and A. Frostegard. 2008. Functional 
robustness and gene pools of a wastewater nitrification reactor: 
comparison of dispersed and intact biofilms when stressed by low 
oxygen and low pH. FEMS. Microbiol. Ecol. 66:167-180. 

12. Mosier, A.C., and C.A. Francis. 2010. Denitrifier abundance and 
activity across the San Francisco Bay estuary. Environ. Microbiol. 
Rep. 2:667-676. 

13. Mou, X., Z. Sun, L. Wang, and C. Wang. 2011. Nitrogen cycle of a 
typical Suaeda salsa marsh ecosystem in the Yellow River estuary. J. 
Environ. Sci. 23:958-967. 

14. Oakley, B.B., C.A. Francis, K.J. Roberts, C.A. Fuchsman, S. 
Srinivasan, and J.T. Staley. 2007. Analysis of nitrite reductase (nirK 
and nirS) genes and cultivation reveal depauperate community of 
denitrifying bacteria in the Black Sea suboxic zone. Environ. 
Microbiol. 9:118-130. 

15. Osaka, T., S. Yoshie, S. Tsuneda, A. Hirata, N. Iwami, and Y. 
Inamori. 2006. Identification of acetate- or methanol-assimilating 
bacteria under nitrate -reducing conditions by stable-isotope probing. 
Microb. Ecol. 52:253-266. 

16. Philippot, L., S. Hallin, and M. Schloter. 2007. Ecology of denitri- 
fying prokaryotes in agricultural soils. Adv. Agron. 96:249-305. 

17. Santoro, A.E., A.B. Boehm, and C.A. Francis. 2006. Denitrifier 
community composition along a nitrate and salinity gradient in a 
coastal aquifer. Appl. Environ. Microbiol. 72:2102-2109. 

18. Schloss, P.D., S.L. Westcott, T. Ryabin, et al 2009. Introducing 
MOTHUR: open-source, platform-independent, community-supported 
software for describing and comparing microbial communities. Appl. 
Environ. Microbiol. 75:7537-7541. 

19. Tamura, K., D. Peterson, N. Peterson, G. Stecher, M. Nei, and S. 
Kumar. 2011. MEGA5: molecular evolutionary genetics analysis 
using maximum likelihood, evolutionary distance, and maximum 
parsimony methods. Mol. Biol. Evol. 28:2731-2739. 

20. Thomas, G., D. Clay, and A. Magurran. 2000. Biodap-ecological, 
diversity and its measurement. Available from http://nhsbig.inhs.uiuc. 
edu/populations/bio-dap.zip. 

21. Yan, T., M.W. Fields, L. Wu, Y. Zu, J.M. Tiedje, and J. Zhou. 2003. 
Molecular diversity and characterization of nitrite reductase gene 
fragments (nirK and nirS) from nitrate- and uranium-contaminated 
groundwater. Environ. Microbiol. 5:13-24. 

22. Yao, S., J. Ni, Q. Chen, and A.G. Borthwick. 2013. Enrichment and 
characterization of a bacteria consortium capable of heterotrophic 
nitrification and aerobic denitrification at low temperature. Bioresour. 
Technol. 127:151-157. 

23. Yoshie, S., N. Noda, S. Tsuneda, A. Hirata, and Y. Inamori. 2004. 
Salinity decreases nitrite reductase gene diversity in denitrifying 
bacteria of wastewater treatment systems. Appl. Environ. Microbiol. 
70:3152-3157. 

24. Zhang, R., V. Thiyagarajan, and P.Y. Qian. 2008. Evaluation of 
terminal-restriction fragment length polymorphism analysis in 
contrasting marine environments. FEMS. Microbiol. Ecol. 65:169- 
178. 

25. Zumft, W.G. 1997. Cell biology and molecular basis of denitrification. 
Microbiol. Mol. Biol. Rev. 61:533-616. 



